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We introduce a new framework of highly-anisotropic hydrodynamics that includes dissipation 
effects. Dissipation is defined by the form of the entropy source that depends on the pressure 
anisotropy and vanishes for the isotropic fluid. With a simple ansatz for the entropy source obeying 
general physical requirements, we are led to a non-linear equation describing the time evolution of the 
anisotropy in purely-longitudinal boost-invariant systems. Matter that is initially highly anisotropic 
approaches naturally the regime of the perfect fluid. Thus, the resulting evolution agrees with the 
expectations about the behavior of matter produced at the early stages of relativistic heavy-ion 
collisions. The equilibration is identified with the processes of entropy production. 
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I. INTRODUCTION 

The experimental results obtained in the heavy-ion ex- 
periments at the Relativistic Heavy- Ion Collider (RHIC), 
in particular the large values of the elliptic flow, are most 
often interpreted as the evidence for a very fast equilibra- 
tion of the produced matter (presumably within a frac- 
tion of 1 fm/c) and for its almost perfect-fluid behavior 

The very fast equilibration is naturally explained 
within a concept that the produced matter is a strongly 
coupled quark-gluon plasma (sQGP) 8]. However, there 
exist other explanations that assume that the plasma is 
weakly interacting. In this case the plasma instabilities 
lead to the fast isotropization of matter, which in turn 
helps to achieve the full equilibration in a short time [9j . 

Recently, several explicit calculations have shown that 
the large values of the elliptic flow and other soft- 
hadronic observables may be successfully reproduced in 
the models that do not assume the very fast equilibra- 
tion. For example, in Refs. (ToJ the stage described by 
the perfect-fluid hydrodynamics was preceded by the free 
streaming of partons (see also [ll|, [l2j), while in Refs. 
[13l [l4j the authors assumed that only transverse de- 
grees of freedom are thermalized (schematic scenarios de- 
scribing the approach towards the full equilibration were 
discussed in this context in Refs. [H, Oil). Such re- 
sults indicate that the assumption of the fast equilibra- 
tion/isotropization might be relaxed. 

It should be also emphasized that the concept of prac- 
tically instantaneous equilibration seems to contradict 



the results of the microscopic models of heavy-ion colli- 
sions. Such models typically use the ideas of color strings 
or color-flux tubes. The system produced by strings is 
highly anisotropic; the pressure in the direction trans- 
verse to the collision axis is usually much larger than the 
longitudinal pressure 1 . 

Similar situation takes place in the Color Glass Con- 
densate (CGC) approach where the distribution func- 
tions are far away from the equilibrium ones. In this 
case the longitudinal momentum distribution is much 
narrower than the transverse one and described by the 
Dirac delta function 6(p\\) at z = [lTlllfl- This approx- 
imation is often used in descriptions of the initial stage 
in nucleus- nucleus collisions (for example see [19|). 

In view of the problems connected with the equili- 
bration and isotropization of the plasma, it is useful to 
develop and analyze the models which can be used to 
describe locally anisotropic systems. In this paper we 
introduce the framework of highly-anisotropic hydrody- 
namics that takes into account dissipation effects. The 
dissipation is defined by the form of the entropy source. 
The latter depends on the pressure anisotropy and van- 
ishes for the isotropic fluid. The proposed model has a 
structure that is very much similar to the perfect-fluid 
hydrodynamics. The main two differences are connected 
with: i) the possibility that the longitudinal and trans- 
verse pressures are different, and ii) the possibility of 
entropy production. By relaxing the assumption about 
the isentropic flow, we generalize our previous formula- 
tions of anisotropic (magneto)hydrodynamics presented 
in[H[2l|]. 

It is important to note that the deviations from equilib- 
rium are naturally described in the framework of viscous 
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1 As usual, the longitudinal direction is defined by the direction of 
the beam. 



(Israel-Stewart) hydrodynamics. However, the region of 
the applicability of viscous hydrodynamics extends to the 
systems that are close to equilibrium. This is reflected in 
the dependence of the transport coefficients on the equi- 
librium variables such as temperature or chemical po- 
tentials. Thus, the viscous corrections are applicable for 
the intermediate, locally almost equilibrated stage (for a 
recent review see [22|). 

In our opinion, the use of viscous hydrodynamics in 
the description of very early stages of the collisions may 
be inadequate — the strong reduction of the initial longi- 
tudinal pressure leads to significant deviations from equi- 
librium (see [23| ). On the other hand, the kinetic models 
that are most suitable for the description of the systems 
out of equilibrium are very much complicated and diffi- 
cult to deal with. Therefore, there is a place for effective 
models which can describe the early non-equilibrium dy- 
namics together with the transition to the perfect-fluid 
regime. Our formulation of the anisotropic hydrodynam- 
ics follows this direction. 

Within our approach, a simple ansatz for the entropy 
source leads to a non-linear equation describing the time 
evolution of the anisotropy in the purely-longitudinal 
boost-invariant systems. The non-linearity implies that 
a possible strong initial anisotropy is eliminated. The 
resulting evolution of the system agrees with the expec- 
tations about the behavior of matter produced at the 
early stages of relativistic heavy-ion collisions. In par- 
ticular, the equilibration of the system is connected with 
the processes of entropy production. 

Although our numerical results are presented for the 
simple one-dimensional system, the proposed formalism 
is general and may be applied to more complicated 2+1 
and 3+1 situations (in a similar way as the perfect-fluid 
hydrodynamics). In addition, different forms of the en- 
tropy source inspired by different microscopic mecha- 
nisms may be analyzed. In our further studies we want 
to explore such rich possibilities. 

Below we assume that particles (partons) are mass- 
less and we use the following definitions for rapidity and 
spacetime rapidity, 
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which come from the standard parameterization of the 
four-momentum and spacetime coordinate of a particle, 

V^ = {E p7 p ± ,P\\) = {p±coshy,p ± ,p x siahy), 
x^ = (t,Xj_, z) = (r cosh 77, x±, r sinhr;) . 

Here the quantity p± is the transverse momentum 



(2) 



Pi. = ypI+Py. 

and r is the (longitudinal) proper time 



(3) 



T = y/t 2 - Z 2 



(4) 

Throughout the paper we use the natural units where 
c = 1 and h = 1. 



II. ANISOTROPIC HYDRODYNAMICS WITH 
DISSIPATION 

A. Energy-momentum tensor and entropy flux 

Our approach is based on the following form of the 
energy-momentum tensor 

T^ = {e + P ± )U"U" -P±g^ - (P± - P^VV , 

(5) 

where e, P±, and P\\ are the energy density, transverse 
pressure, and longitudinal pressure, respectively. In the 
special case of the isotropic fluid, where P± = A = P, 
we recover the form of the energy-momentum tensor of 
the perfect-fluid hydrodynamics. The four-vector U^ de- 
scribes the hydrodynamic flow, 



EP = 7(l,v), 7 =(l-« 2 ) 
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while V^ defines the direction of the longitudinal axis 
that plays a special role due to the initial geometry of 
the collision. The four-vectors U^ and V^ satisfy the 
following normalization conditions 



U- 



1. 



V- 



1, 



U-V = 0. 



(7) 



In the local-rest- frame (LRF) of the fluid element we have 
U^ = (1,0,0,0). In the same reference frame the four- 
vector V^ that satisfies Eqs. and defines the longitu- 
dinal direction has the form V* = (0, 0, 0, 1). The forms 
of C/ M and V** in other reference systems are obtained by 
the Lorentz boosts. 

In LRF the energy-momentum tensor has the diagonal 
structure, 
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Hence, as expected, the formula (|S|) allows for different 
pressures in the longitudinal and transverse directions. 

In addition to the energy-momentum tensor ([5]) we in- 
troduce the entropy flux 



aU", 



(9) 



where a is the entropy density. We assume that e and 
a are functions of P± and P« . In particular, since we 
consider massless partons here, the condition T M = 
gives 



e = 2P, 



P 



(10) 



We note that the form of the energy-momentum tensor 
([5]) resembles the form used in relativistic magnetohydro- 
dynamics [2J, (25| . In that case the anisotropy is induced 
by the presence of the magnetic field. At the early stages 
of heavy-ion collisions we have similar situation — there 



exist strong color magnetic and electric longitudinal fields 
(Glasma following CGC [26]) which polarize the medium. 
The explicit inclusion of the fields should be one of the 
first tasks connected with the generalization of the pre- 
sented framework. The first steps in this direction but 
for the isentropic case were made in |2l| . 



B. Evolution equations 

The dynamics of the system is governed by the equa- 
tions expressing the energy-momentum conservation and 
the entropy growth (the second law of thermodynamics), 



d„T^ = 0, 
d^ = £. 



(11) 
(12) 



Here the function £ describes the entropy source. The 
form of £ must be treated as the assumption defining the 
dynamics of the anisotropic fluid. Besides the condition 
£ > it is natural to assume that £ = for P± = P\\ . In 
this way, in the case where the two pressures are equal, we 
recover the structure of the perfect-fluid hydrodynamics. 
We note that £ is the internal source of the entropy, i.e., 
£ describes the entropy growth due to the equilibration 
of pressures in the system. 

In the following we shall treat £ as a function of P± 
and P|| . In this way, Eqs. (fTTj) and (fT2j) form a closed 
system of 5 equations for 5 unknown functions: three 
components of the fluid velocity, P±, and Pm. The pro- 
jections of Eq. (fTTj) on U u and V„ gives 
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C. Anisotropic momentum distribution 
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In our previous paper [2l| we showed that the structure 
(O-© follows from the following form of the distribution 
function 
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Here the two parameters Aj_ and A|| may be interpreted as 
the transverse and longitudinal temperature. The form 
(fT5"|) is valid in the local rest frame of the fluid where 
W = (1,0,0,0) and V = (0,0,0,1). The explicitly 
covariant form of the distribution function (|15[) has the 
form 



!=f f ^-ur-^vr ^\ (16) 



In this paper we consider the exponential distribution 
function which in the local rest frame has the form 



/ = 3o exp - 
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Equation p7|) may be regarded as the generaliza- 
tion of the Boltzmann equilibrium distribution where 
Aj_ = A|| = T. The parameter go is the degeneracy factor 
connected with internal quantum numbers. Having in 
mind the fact that the initially produced matter consists 
mainly of gluons we obtain 



.9o = 16. 



(18) 



Using the covariant form of (|17[) in the definition of 
the energy-momentum tensor, 
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and in the definition of the entropy flux 2 



d 3 p p^ 



I 



1 -In 



(19) 



(20) 



we obtain Eqs. ((5} and ((9]). 

The energy density, transverse pressure, longitudinal 
pressure, and entropy density are obtained from the fol- 
lowing integrals 
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D. Pressure anisotropy 

Equations ([2"Tj) - (T2"4")) allow us to express all thermody- 
namic quantities 3 in terms of Aj_ and An. Thus, instead 



2 The formula H20H assumes the classical Boltzmann statistics. It 
may be generalized to the case of bosons or fermions in the stan- 
dard way. 

3 We continue to call 6, P± , Pu , and a the thermodynamic quan- 
tities although, strictly speaking, these quantities do not refer to 
the equilibrium state. Similarly, we call Ax and Am the transverse 
and longitudinal temperature. The reason for this terminology 
is close similarity to the equilibrium variables. 



2(1 



10 



L 



R(x) 



R'(x) 



R(l) 



3 go 



R' (1)=0 



0.0 0.5 1.0 1.5 2.0 2.5 



10 r 



11.(1 



P„ R (x) I 3 - 2 x R' 
P ± R (x) / 3 + x R' (x) 



(1.5 



1.0 



1.5 
X 



2.0 



2.5 



3.0 



FIG. 1: (Color online) The ^-dependence of the function R(x) 
and its derivative R' (x). 



of Pj_ and P\\ we may switch to A_l and A|| . It turns out, 
however, that the most useful two independent parame- 
ters are the entropy density a and the variable x defined 
by the expression 



(25) 



Treating a and x as the two independent thermodynamic 
variables (instead of P± and Pn or instead of Aj_ and Ai|) 
we obtain the compact expressions: 
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where the function P(x) is defined by the formula [2J| 4 
3s , oa;~ t 
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The symbol R'(x) denotes the derivative of R(x) with 
respect to x. For x = 1 we find P'(x) = and, as 
expected, P± = P\\ . The ^-dependence of the function 
R(x) and its derivative R'(x) is shown in Fig. [TJ 

In Fig. [5] we show the ratio of the longitudinal and 
transverse pressure plotted as a function of x (solid 



4 Note that for x < 1 the function (arctan y/x — l)/y/x — 1 should 
be replaced by (arctanh\/l — x)/ \/l — x 



FIG. 2: (Color online) The ^-dependence of the ratio of the 
longitudinal and transverse pressure (thick solid line) and its 
approximation with the function x~ 3 ' 4 (thick dashed line). 



thick line). The P||/Pj_ ratio is determined by the x- 
dependence of the function R(x) and its derivative. To a 
good approximation one finds 



5l 
Pi 



-3/4 



(30) 



Thus x may be treated as the direct measure of the pres- 
sure anisotropy. We note that realistic initial conditions 
in heavy- ion collisions give Pii/Pj_ <C 1, which corre- 
sponds to x >■ 1. 



III. PURELY-LONGITUDINAL 
BOOST-INVARIANT MOTION 

A. Implementation of boost-invariance 

For purely longitudinal and boost-invariant motion we 
may write 



and 



P M = (cosh??, 0,0,sinh?7), 



V^ = (sinh?7, 0,0,cosh7y). 
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The boost-invariant character of Eq. (|TB|) is immediately 
seen if we write the explicit expression for p ■ U and p ■ V, 
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We also obtain 
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(34) 



which leads to the equations 



we^ir = o, 

tV^V = U v , 
0„V" = 0. 



(35) 



We note that the boost-invariance requires that all scalar 
quantities such as s, pr, or Pl do not depend on space- 
time rapidity rj. 



B. Boost-invariant equations of motion 

In the considered case, the energy-momentum conser- 
vation law (|13p is reduced to the equation 



de __ e + Pj| 
dr t 

while the entropy conservation yields 

da 1 _ £ 

adr t a 



(36) 



(37) 



Equation (TT41 is automatically fulfilled if the thermody- 
namic variables do not depend on r\. 

By changing to the x variable we may rewrite Eq. (|36[) 
in the form 



R\x) 



dx 
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Before we proceed further with the analysis of the dis- 
sipative flow where £ > 0, it is useful to consider the 
non-dissipative flow where £ = 0. In this case, from Eq. 
(|37p we recover the Bjorken solution a = ctoTo/t and the 
right-hand-side of Eq. (j3"5)l vanishes. This implies that 
either x = 1 (in which case R'(l) — 0) or x — xqt 2 /tq 
(in which case dx/dr — 2x/t). The parameters ctq, to, 
and xq are arbitrary constants here. The case x = 1 cor- 
responds to the standard, perfect-fluid hydrodynamics. 
The case where x = xqt 2 /t§ was examined in [20L l2llj . 



Ansatz for E 



Equations (j3"7) and ([55} may be solved only if the de- 
pendence of the function £ on the variables a and x is 
given. The functional form £(<r, x) must be delivered as 
the external input for our calculations. 

The simplest ansatz for £ that has the correct dimen- 
sion and satisfies the general conditions that £ > and 
£(<7, x = 1) = has the form 
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Here r eq is a timescale parameter 



nuit / eq io ct u-ii~u.c<Qv-.ciic< \jcAii.cxiL±L\^ii^i. . ±. ne natural ieaxure 
of (131))) is the fact that £ is proportional to a, hence (|3"9")l 
does not destroy the scale invariance of the perfect-fluid 



hydrodynamics, that allows for multiplication of a in the 
evolution equations by an arbitrary constant. Moreover, 
£ defined by (|39|) stays constant if Aj_ and A|| are inter- 
changed. 

Substituting (|39|) in (f38|) leads to the ordinary differ- 
ential equation for x only, 



dx 

where we have defined 
H(x) = 
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D. Results 

We solve Eq. (|40|) numerically with the initial con- 
dition x = xq set at r = ro = 0.2 fm. The results 
of the microscopic models suggest that Pm <C P± at 
the early stages of the collisions, thus we first con- 
sider the case Xq = 100. For completeness, in the end 
of this Section we also show the results obtained with 
xq = 0.01. The time evolution is studied in the time in- 
terval: 0.2 fm < r < 10 fm. 

In Fig. |3jwe show the time dependence of various phys- 
ical quantities obtained with xq = 100 for three differ- 
ent choices of the relaxation time: r cq = 0.25 fm (solid 
line), r cq = 0.5 fm (dashed line), and r cq = 1.0 fm (dot- 
ted line). 

Figure [3] (a) shows the time dependence of the asym- 
metry parameter x. We observe a fast changes of x in the 
initial stages of the evolution. Such changes, depending 
on T Gq , are caused mainly by the fact that H (x) behaves 
like 6a; 3 / 2 for large values of x. Hence, large initial values 
of x imply also large (but negative) values of the deriva- 
tive dx/dr at t = To . We discuss this behavior in greater 
detail in Sect. IVAl 

The behavior shown in Fig. [3] (a) indicates also that 
x ps 1 for t > 2r cq . This is a desired effect showing that 
the system approaches the local equilibrium state. The 
way how x approaches unity is described in more detail 
in Sect. IV Bl where the approximate analytic solution is 
presented. 

In Fig. [3] (b) we compare the time evolution of the 
entropy density obtained from Eq. (|38|l with the Bjorken 
solution 
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Here cro is the initial value of the entropy density. We 
note that the specific value of <7o is irrelevant for our 
analysis, since the entropy equation is invariant with re- 
spect to the multiplication of a by an arbitrary constant. 
The amount of the entropy produced in the 
regime described by the anisotropic hydrodynamics 
depends in our case on the relaxation time. For 
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FIG. 3: (Color online) (a) The time dependence of the asym- 
metry parameter x for three different choices of the relaxation 
time: r cq = 0.25 fm (solid line), r cq = 0.5 fm (dashed line), 
and r cq = 1.0 fm (dotted line), (b) Entropy density divided 
by the corresponding values obtained in the Bjorken model, 
(c) Ratio of the longitudinal and transverse pressures shown 
as a function of the proper time. All results are obtained with 
the initial asymmetry xq = 100. 



7~cq = 0.25, 0.5, and 1.0 fm the entropy increases by 
about 70%, 85% and 105%, respectively. For r > r cq 
the ratio (err) / (<jqTo) saturates indicating that the flow 
attains the form of the Bjorken flow. This behavior shows 



FIG. 4: (Color online) The same as Fig. [3] but for the initial 
condition x(jo) = xo — 0.01. 



again that our framework may be used to model the tran- 
sition between the highly anisotropic initial behavior and 
the perfect-fluid stage. 

In Fig. [3] (c) we show the time dependence of the ra- 
tio of the longitudinal and transverse pressure. Again, 
we show three different time evolutions corresponding to 
three different relaxation times. For r 3> r cq the ratio 
approaches unity and the two pressures become equal. 

Figure 0] shows the same time evolutions as Fig. [3] 



but with the initial condition xq — 0.01. Our main re- 
mark here is that the dynamics of the system governed 
by Eqs. (|38]l and (|37|) again leads to the equilibration of 
the system. The entropy production in the anisotropic 
phase is similar to the previous case. Interestingly, the 
time dependence of a: and P\\/P± is not monotonic in this 
case, but anyway P± ~ Pu for sufficiently large evolution 
times. 



IV. CONCLUSIONS 

In this paper we have introduced the new framework 
of highly- anisotropic hydrodynamics with strong dissipa- 
tion. The effects of the dissipation are introduced by the 
special form of the internal entropy source. The source 
depends on the pressure anisotropy and vanishes for the 
isotropic systems to guarantee that the perfect-fluid be- 
havior is reproduced for the locally equilibrated system. 

With a simple ansatz for the entropy source satis- 
fying general physical requirements, we have obtained 
a non-linear equation describing the time evolution of 
the anisotropy parameter x. The non-linearity causes 
that the initial large (or small) anisotropy parameter ap- 
proaches asymptotically unity. The rate at which the 
equilibrium is reached depends on the relaxation-time pa- 
rameter T eq . We think that with a suitable chosen value 
of T cq , our approach may be useful to model the fast equi- 
libration of matter expected in heavy-ion collisions. In 
particular, it offers an attractive option for modeling the 
continuous equilibration of pressures. 

The dynamics of anisotropic fluid determines the 
changes of the microscopic distribution function f(x,p). 
Thus, various calculations done so far for the systems 
in equilibrium may be repeated for the non-equilibrium 
case. In this way the effects of the non-isotropic dynam- 
ics may be analyzed and, more importantly, verified in a 
very straightforward way. 

Our numerical results have been presented for a sim- 
ple one-dimensional system. Nevertheless, the proposed 
formalism is general and may be applied to more com- 
plicated 2+1 and 3+1 situations. In addition, different 
forms of the entropy source inspired by different micro- 
scopic mechanisms may be analyzed. In this context it 
is interesting to search for the hints coming from the 
ADS/CFT correspondence. In our further studies we 
want to explore such rich possibilities. 
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FIG. 5: (Color online) Function H(x) and its approximations 
in the regions x ^> 1 (a) and \x — 1| <C 1 (b) . 



see Fig. [5] (a). This leads to the equation 
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which has the following analytic solution 
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V. APPENDICES 

A. Analytic solutions for x S> 1 

For very large arguments H (x) may be approximated 
by the formula, 



H(x)^6x 3/2 , 



(43) 



We note that large initial values of x imply very large 
negative values of the derivative dx/dr, since directly 
from Eq. (14"4"|) one finds 



dx 
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B. Analytic solutions for \x — 1| -C 1 

In the region \x — 1| <C 1 we may use the following 
approximation, see Fig. [5] (b) , 
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H(x)*-(x-l) + —(x-l) 2 



(47) 



If x is close to 1 we take the first term in the series and 
obtain 



dx 2x 



15 



(3-1). 



I cq 



The solution of this equation has the form 



x = —5- exp 



15t 



.4 



225, 



-Ei 



15r 

4r 0Q 



(48) 



15t 

4^T 



CC| 

(49) 



where Ei(ir) is the exponential integral function and A is 
an arbitrary integration constant. Using the asymptotic 
expansion of Ei(x) for x ^> 1, 



Ei 



. /15t 



eq 



exp 




. (50) 



one obtains 



lim x(t) = 1. 

t>-r oq 



(51) 



Similar, but much more involved calculations may be 
done for the case where the second term in the series 
(1471) is included. 
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